Oscillatory dynamics and organization of the vortex solid in YBa2Cu3C>7 single crystals 
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We report on the degree of order of the vortex solid in YBa2Cu307 single crystals observed in ac 
susceptibility measurements. We show that when vortices are "shaken" by a temporarily symmetric 
ac field they are driven into an easy-to-move, ordered structure but, on the contrary, when the 
ac field is temporarily asymmetric, they are driven into a more pinned disordered state. This is 
characteristic of tearing of the vortex lattice and shows that ordering due to symmetric ac fields is 
essentially different from an equilibration process or a dynamical crystallization that is expected to 
occur at high driving currents. 



PACS numbers: 74.60. Ge, 74.60. Jg 

Vortex lattices (VL) in type II superconductors provide 
a unique system with tuneable parameters for exploring 
driven dynamic phases. A characteristic feature of driven 
dynamics in quenched disorder is the existence of a de- 
pinning transition. Experimental |l| and theoretical ||] 
results have shown evidence of a two-step process in the 
depinning transition of the VL as the driving force is in- 
creased: first it undergoes plastic flow, where neighboring 
parts of it move at different velocities, and then, above 
a threshold force Ft, a dynamic crystallization occurs as 
proposed by Koshelev and Vinokur Q (K-V) where all 
the vortices move at the same average velocity j| . With 
the use of fast current ramps and ac techniques Q, the 
attention has been focused on the evolution of the VL 
as it dcpins and starts moving. Though the idea of this 
annealing process is accepted in steady driven VL's, a 
number of new puzzling phenomena were observed with 
these techniques that cannot be ascribed to it such as 
complex history dependent dynamic response and mem- 
ory effects 0-0]- To explain these phenomena, it has 
been proposedlq] that the overall vortex structure surges 
as a balance between the injection of disordered vortices 
through surface barriers and an annealing by the driv- 
ing current. However, this seems not to be the case in 
the high T c twinned YBa2Cu307 (YBCO) single crystals, 
where similar effects have been observed, but only bulk 
pinning appears to be relevant H . In the latter investiga- 
tion, it was found that an "easy- to-move" more ordered 
VL can be stabilized by the shaking movement induced 
by a temporarily symmetric ac field (e.g. sinusoidal). 
Very recently, it has been suggested that this reorder- 
ing may be a consequence of a current-assisted transition 
from a supercooled metastable phase to the stable equi- 
librium phase |@,@]. The fluctuation energy produced in 
the sample under cyclic field variation would be a key 
element within this scenario . 

In this Letter, we show that if the shaking field is tem- 
porarily asymmetric (e.g. sawtooth) a disordered VL 
can be recovered. The "hard-to-move" state so obtained 
would have a high density of topological defects due to 
the plastic flow of vortex bundles sliding one past an- 



other. The fact that the degree of order is strongly de- 
pendent on the temporal symmetry of the ac field would 
indicate that the reordering observed in Ref. f| is due 
neither to a K-V transition nor to a current-assisted equi- 
libration but to a process of a fundamentally different 
kind which involves repeated interaction of neighboring 
vortices. 

We also propose that when the VL is pushed in and 
out of the sample with an asymmetric ac field, it is torn 
up in a "ratchet-like" fashion as a consequence of the 
different Lorentz forces involved in the ramp up and ramp 
down of the field. Important related results show that 
these plastic distortions are partly reversible when the 
temporal evolution of the ac field is reversed. 

The results shown here are for a single crystal sample 
of YBCO H (dimensions 0.56 x 0.6 x 0.02 mm 3 ) with T c 
= 92 K and ST C = 0.3 K determined by ac susceptibility 
(h ac =l Oe) at zero dc field. The same sample was used 
in previous studies on plasticity and memory effects M. 
Measurements were done with the standard mutual in- 
ductance technique, with the ac field parallel to the c axis 
of the sample. The in-phase and out-of-phase (90°) sig- 
nals were collected with a lock-in amplifier. Data were 
recorded for an angle 9 = 20° between the dc applied 
field, Hdc, and the c axis to avoid the Bose-glass phase 
f|. The dc field was oriented out of all twin boundaries 
simultaneously (see inset in Fig. 1). 

Fig. 1 presents ac susceptibility measurements that 
correspond to different thermomagnetic histories of the 
sample. Measurements were performed by varying T (in 
the direction indicated by the arrows) at a rate of 0.2 
K/min, at fixed H^ c (3 kOe), h ac (1.6 Oe) and frequency 
(10.22 kHz) (default parameters). Dashed curves were 
obtained in the usual field cooled procedure (F ac C C) 
while solid ones were obtained on warming after zero (ac 
and dc) field cooling (ZFC W) (the ac field is turned on 
after the dc field has reached its final value Jl(J). We 
believe that the ZFC W case corresponds to a disordered 
state because of the plastic motion of vortices as they 
penetrate into the sample ||. In addition, in the F ac C C 
case, the VL is partially annealed by the shaking move- 
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FIG. 1. x"(T) (a) and X '(T) (b) for the YBCO single crys- 
tal with different thermomagnetic histories. Arrows indicate 
the direction of the temperature sweep. Dashed curves were 
obtained on cooling, F ac C C. Solid curves were obtained on 
warming after cooling with neither ac nor dc field applied, 
ZFC W. Application of 10 cycles of a triangular ac field (6.5 
Oe, 10 kHz) to a F QC C C state at 85. 2K (A) leads to B and 
upper curves (open circles). Application of 10 5 cycles of a 
sawtooth ac field (6.5 Oe, 10 kHz) to a F ac C C state at 85. 2K 
(A) leads to C and lower curves (solid circles). In all cases 
hac | c — 1.6 Oe and Hdc = 3 kOe is applied at 6 = 20° from 
the c axis (see lower inset in the upper panel.) 



merit induced by the ac field during the cooling process 
Jil"| ]. As discussed previously ||, the history effects are 
observed because of a change in the mobility of the VL 
and the fact that the inner disordered state in the ZFC W 
case can sustain a higher current without vortex move- 
ment thereby enhancing the shielding \x'\ and reducing 
the dissipation x" (for small measuring ac fields). 

Through the F ac C C procedure one can obtain the 
most ordered state attainable with the measuring ac field. 
However, a more ordered VL can be obtained if the VL 
is shaken more intensely with a higher ac field. To illus- 
trate this, we prepared a F ac C C state at 85.2 K (point 
A) and applied for a short time (ten seconds) a triangular 
ac field with h ac — 6.5 Oe and /= 10 kHz. We turned 
the 6.5 Oe ac field off and measured the ac susceptibility 
with our default parameters on warming. The ordering 
of the VL is inferred from the unambiguous changes of 
x' and x" (point B and upper curves in Fig. 1, open 
circles). 

Our first important result surges when contrasting 
these measurements with the ones obtained after the ap- 



plication of an asymmetric ac field under the same con- 
ditions. We prepared again a F ac C C state at 85.2 K 
(point A) and applied during ten seconds a sawtooth ac 
magnetic field (h ac = 6.5 Oe, / = 10 kHz). The mag- 
netic field waveform has a rising time, i+, and a shorter 
falling time, t_, (see inset in Fig. 1(a)) so that a given 
current forces vortices, for example, towards the center of 
the sample during the longer time interval, but a higher 
current forces them out during the rest of the cycle. We 
turned the 6.5 Oe ac field off and then measured the ac 
susceptibility. The result is indicated as point C in the 
lower curves of Fig. 1(a) and (b), where it is clearly seen 
that \x'\ increases and x" decreases. On further warming, 
we found that the measured susceptibility (solid circles) 
goes close to the disordered ZFC W curve. The difference 
in the response after the application of a symmetric or an 
asymmetric ac field is striking. The sawtooth ac field has 
brought the VL to a more pinned disordered state. It is 
worth noting here that one can go from point C to point 
B, or vice versa, by properly applying a high amplitude 
symmetric or asymmetric ac field as explained above. 

Fig. 2 exhibits clear evidence that both the order- 
ing and disordering of the VL sketched in Fig. 1 are 
number-of-cycles dependent. The same information can 
be extracted from x' an d x" ■ As the relative change in 
x" is larger than in x'i we only present x" data. The 
initial states are obtained after exposing the sample to 
10 4 cycles of sawtooth ac field (initial disordered state in 
Fig 2(a)), plus 10 5 of triangular ac field (initial ordered 
state in Fig. 2(b)) g§. 

Fig. 2(a) shows the dynamics of the VL's ordering. 
Starting from a disordered hard-to-move state, the sam- 
ple is cycled with a symmetrical (sinusoidal, triangular 
or square) ac field (h ac = 6.5 Oe, f= 10 kHz). After N s 
cycles, the ac field is turned off and the state of the VL is 
inspected by measuring Xac with a smaller probe ac field 
(default parameters above). Before repeating the proce- 
dure for a new value of N s , the VL is brought again to the 
initial disordered state. It is seen that the dependence of 
Xac on N s is approximately logarithmic, and saturation, 
if any, occurs beyond 10 6 cycles. 

The evolution of the disordering was inspected in a 
similar experiment at the same temperature. In Fig. 
2(b) we show Xac as a function of N a asymmetric ac 
field cycles applied to an initially easy-to-move ordered 
lattice. The asymmetry of the waveform is defined as 
a = (i + /<_). In all cases, the preparation of the ini- 
tial ordered state before the application of the next N a 
asymmetric ac field cycles erases any eventual effects of 
the small ac probe used in the Xac measurement. 

Ordering due to the application of a driven current or 
an oscillatory perturbation has been frequently observed 
and proposed to explain some unusual behaviors 
Some of these results have been interpreted as an equili- 
bration process |7],|| or as a combined effect of the injec- 
tion of disordered vortices through the sample edges and 
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FIG. 2. Oscillatory organization of vortices, (a) Ordering 
of an initially disordered VL: x vs - N a for triangular (solid 
triangles) , sinusoidal (open circles) and square (solid squares) 
ac field symmetric waveforms with h ac = 6.5 Oe and / = 10 
kHz. (b) Disordering of an initially ordered VL: x" vs - N a for 
different asymmetries a of the field (see text), h ac = 6.5 Oe 
and / = 10 kHz. a = 1 corresponds to a triangular waveform 
and ol = oo to a sawtooth waveform. Measuring ac field: 1.6 
Oe, 10.22 kHz. Lines are guides to the eye. 

a crystallization assisted by the driving current ||. In 
the former case, it has been proposed |^ that metastable 
to stable transformations occur in the regions of the sam- 
ple where the local energy dissipation exceeds a threshold 
value. 

In our experiments, it is clear that an asymmetric ac 
field will also produce an oscillatory perturbation; but 
then it is legitimate to inquire into why the VL disor- 
ders. The discrepancy that this implies suggests that the 
ordering process is unlikely to be a K-V crystallization 
and that it is surely not an equilibration transition where 
energy dissipation in the sample will always tend to order 
the VL g. 

In the K-V scenario Q , ordering is only achieved when 
the Lorentz force exceeds a certain threshold value, Ft, 
while a plastic disordered motion occurs for forces below 
it. The induced current density in our experiments de- 
pends directly on the sweep rate of the applied ac field 
which, on the other hand, is proportional to its frequency. 
Dynamical reordering is observable, at least, for frequen- 
cies in the range of 0.1 Hz to 300 kHz for sinusoidal, trian- 
gular and even for the high sweep rate square waveform. 
Disordering shown in Fig. 2(b) is clearly discernible for a 
frequency (10 kHz) well within the range of frequencies in 
which ordering is achieved and for asymmetries as small 
as 3. This observation would indicate that for equal val- 
ues of current density circulating in the sample one can 
either order or disorder the VL depending only on the 
asymmetry of the ac field and that, remarkably, ordering 
can be achieved with lower current densities (e.g., at 0.1 
Hz) than those that are able to disorder the VL as shown 
in Fig. 2(b). On account of this, the K-V transition can 
be discarded as a possible explanation of our observa- 



tions. In addition, the dynamically generated disorder 
with the asymmetric ac field is distinctive of tearing of 
the VL, which suggests that the driving force is below 
Ft and that the flow of vortices is plastic (first step in 
the depinning process In these conditions, there 

are regions of the VL where the strain is large enough to 
cause phase slips and the number of vortices that partici- 
pate in the movement strongly depends on the magnitude 
of the applied force, ie. on the induced current. 

If the ac field is asymmetric, the vortices that move 
when the field is ramped up will not be the same ones 
that move when the field is ramped down. Then, after 
the application of one cycle of an asymmetric ac field, 
there should be a net displacement between vicinal bun- 
dles of vortices. This type of movement will cause lattice 
distortions that destroys the long range order of the VL 
and leads to the increase of \x'\ and the decrease of \" 
for small ac fields. Note that these distortions are cu- 
mulative and increase with the number of cycles leading 
to a rapid saturation of Xac after ~ 300 cycles for the 
sawtooth waveform. When a is decreased, the disorder 
also decreases and the number of cycles needed to reach 
maximum disorder increases. 

Recovering an ordered VL requires a larger number 
of cycles of a symmetric ac field of the same amplitude. 
As seen in Fig. 2(a), after 10 6 cycles \' an d x" still 
vary approximately as In AT. It is the overall symmetry 
of the ac field that drives vortices back and forth, with 
no net displacement between them. The effect of this 
shaking movement of vortices is to locally order the VL. 
Apparently, repeated interactions of vortices with their 
nearest neighbors seems to enhance ordering and favor 
the healing of defects. 

A curious behavior is observed when a is reduced: the 
change in x" is non-monotonous in the number of cycles. 
This is more clearly seen for the smaller asymmetries. 
Initially, x" rapidly decreases but after ~ 10 4 cycles this 
tendency is reversed and x" slowly starts to increase. 
This phenomenon could be explained through a process 
involving dynamic and static friction. The argument 
could be as follows. At the beginning, the VL disorders 
as described above as a consequence of the asymmetry of 
the ac field. However, as a is small, the bundles of vor- 
tices that just start moving when subjected to the larger 
Lorentz force (higher field sweep rate), were at the static 
limit in the preceding part of the field cycle (lower field 
sweep rate). Once they move, it is established a new 
state determined by dynamic friction (which presumably 
is smaller than the static one) in which these bundles of 
vortices move coherently with the rest. If this picture 
is correct, it would take about ~ 10 4 cycles to reach a 
steady state in which vortices move coherently and the 
repeated interaction between them, that allows ordering, 
dominates in what Xac is concerned. 

Finally, we address the question whether the distor- 
tions produced by the asymmetric ac field can be released 



3 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



r= 85.2 K 

H =3 kOe; h - 1.6 Oe 




10" 



FIG. 3. Releasing distortions in the VL. x" as a function 
of the number N+ of sawtooth field cycles applied to an ini- 
tially ordered VL. After applying N+max cycles, the polarity 
is inverted as shown in the inset. Nt is the total number of 
cycles N+ + N-. Curves are vertically displaced for clarity. 
Measuring ac field: 1.6 Oe, 10.22 kHz. Lines are guides to 
the eye. 



if we make a temporal inversion of the ac field. If vortex 
movement were completely reversible and we applied N+ 
cycles of an asymmetric ac field, we would recover the 
initial configuration of the VL after applying the same 
number of cycles of an asymmetric ac field with the po- 
larity inverted. The results in Fig. 3 suggest that a more 
ordered configuration can be recovered whenever the dis- 
tortions generated are below certain level. As shown in 
the inset of Fig. 3, we first apply a given number of 
cycles, N + , of sawtooth ac field to an initially ordered 
VL. As in Fig. 2, \" decreases as N + increases. After 
N+ m ax cycles, we invert the polarity of the ac field so 
that the next cycles represent a "temporal inversion" of 
the previously applied ac field. At this point it is ob- 
served that x" starts increasing and, when the number 
of cycles of the inverted waveform, N-, exceeds N +max , 
x" decreases again. We have followed this procedure for 
N +max = 5, 10, 20, 40 and 256. When N +max is high 
enough to highly distort the VL and saturation of x" 
is reached (N +max = 256) almost no reversibility is ob- 
served. These results further support the scenario pro- 
posed of ratchet-like tearing of the VL by a temporarily 
asymmetric ac field. 

To conclude, we have presented susceptibility measure- 
ments in a twinned YBCO single crystal after shaking the 
VL with an ac field. We find a dramatic change of the 
VL response when the asymmetry of the shaking ac field 
is increased. While a symmetric ac field orders the VL 
an asymmetric ac field produces the opposite effect. The 
experiments described demonstrate that ordering is not 
due to a K-V transition or to an equilibration process. 
The ordering of the VL appears to be a consequence of 



the repeated interaction of neighboring vortices while the 
disorder produced by an asymmetric ac field would be 
related to tearing caused by the different Lorentz forces 
involved in the ramp up and down of the ac field. We 
hope these results will stimulate both experimental and 
theoretical work on the dynamics of VL's under an os- 
cillatory perturbation to further test the hypotheses pre- 
sented here. 
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